[W,size_t] = stft_2(x,fs,4); %takes stft
W_abs = 10*log10(abs(W));
C(:,:,i) = chromatic(x',y, 1024); %takes chroma/NPCP
c_sum(:,i) = sum(C(:,:,i), 2);
%generate plots for IF
h = figure;
xticklabels = 0:60:ceil(length(x)/y);
yticklabels = 0:500:ceil(y/2);
xticks = linspace(1, nFrames, numel(xticklabels));
yticks = linspace(1, frameSize, numel(yticklabels));
axis xy;
colormap jet
imagesc(iFz)
caxis ([1200 3200])
c = colorbar;
c.TicksMode = 'manual';
c.Ticks = [-2000 0 2000 4000 6000 8000];
c.TickLabelsMode = 'manual';
c.TickLabels = {'-2000', '0', '2000', '4000', '6000', '8000'};
set(gca, 'YDir', 'normal');
xlabel('Time (seconds)');
ylabel('Frequency (Hz)');
title({'Instantaneous frequency of track: ' filename(i).name});
saveas(gca,['Chroma(y)', filename(i).name , '.png']);
close(h);
%generate chroma plots
h=figure;
C_abs = 10*log10(abs(C));
imagesc(flipud(C_abs));
ax=gca;
ax.YTickLabel=fliplr({'A ','A#','B ','C ','C#','D ','D#','E ',...
'F ','F#','G ','G#'});
ax.YTick=linspace(1,12,12);
colormap 'jet';
title({'Chroma :'; filename(i).name});
xlabel('frames');
ylabel('Note');
colorbar;
caxis([-60 0])
saveas(gca,['NPCP' filename(i).name '.png']);
close(h);
end
% [x1,fs] = audioread('1.wav');
% x2 = audioread('3.wav');
%
% [C1] = chromatic(x1',fs,1024);
% [C2] = chromatic(x2',fs,1024);